Valley separation in graphene by polarized light 
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We show that the optical excitation of graphene with polarized light leads to the pure valley 
current where carriers in the valleys counterflow. The current in each valley originates from asym- 
metry of optical transitions and electron scattering by impurities owing to the warping of electron 
energy spectrum. The valley current has strong polarization dependence, its direction is opposite 
for normally incident beams of orthogonal linear polarizations. In undoped graphene on a substrate 
with high susceptibility, electron-electron scattering leads to an additional contribution to the valley 
current that can dominate. 
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Graphene, one-atom-thick layer of carbon with the 
honeycomb crystal lattice, has been attracting rapidly 
growing attention due to its unique electronic properties. 
Zero band gap and zero effective electron and hole masses 
as well as high enough mobility make it perspective for 
fundamental and applied research!^! The electron exci- 
tations in graphene are similar to massless Dirac fermions 
with the cone points situated at the points K and K' 
of the Brillouin zone. The interplay of two equivalent 
valleys gives rise to new transport and optical phenom- 
ena, which are absent in systems with simple electron 
dispersion, an d un derlies the novel research field called 
"valleytronics" PE1 In multivalley structures, one can in- 
dependently control the carriers in different valleys and 
construct peculiar electron distribution where particles 
in the valleys flow predominantly in different directions.^ 

Previous research of valley-dependent transport in 
graphene was focused on the manipulation of charge 
carriers by static electric field. It was demonstrated 
that the electric field may induce valley-polarized cur- 
rent in a graphene point contact with zigzag edges, 5 
graphene layer with broken inversion symmetryJ-3 bilayer 
graphene,^ or if the structure is additionally illuminated 
by circularly polarized radiation!^ It was also proposed 
in Ref. [10] that valley currents can be induced in meso- 
scopic graphene rings by asymmetrical monocycle elec- 
tromagnetic pulses. Here, we show that the valley sepa- 
ration can be achieved in a homogeneous graphene layer 
by pure optical means. We demonstrate that the inter- 
band excitation of graphene by linearly polarized light 
leads to the electron current in each valley, which direc- 
tion is determined by the light polarization. The partial 
photocurrents (y = ± for the valleys K and K 1 ', re- 
spectively) in the ideal honeycomb structure are directed 
oppositely, so that the total electric current jW + 
vanishes. We also briefly discuss optical and transport 
methods to reveal the pure valley current. 

Phenomenologically, the emergence of the valley pho- 
tocurrent is related to the low point-group symmetry of 
individual valleys. Despite the fact that the crystal lat- 



tice of flat graphene is centrosymmetric, the valleys K 
and K' are described by the wave vector group D$h lack- 
ing the space inversion, see Fig. 1. The group D^h allows 
for the photocurrent induced by normally-incident lin- 
early polarized light. The polarization dependences of 
the current components in the valley K are given by 

= X(e 2 x - e 2 y )I , jW=-2 X e x e y I. (1) 

Here, ^ is a parameter, e x and e y are components of 
the (real) light polarization unit vector e, / is the in- 
tensity of incident light, and the x axis is chosen along 
the T-K line, Fig. la. The photocurrent in the valley K' 
is obtained from Eq. (fTl) by the replacement x — > —x, 
which gives j^' = — We note that the absence of 
a total electric current at normal incidence of radiation 
is in agreement with the symmetry arguments allowing 
for a photocurrent in noncentrosymmetric systems only. 
At oblique incidence of the radiation, a net cur rent i n 
graphene may arise due to the photon drag effect lilllU 

The microscopic model of pure valley current genera- 
tion is based on the trigonal warping of energy spectrum 




FIG. 1: (color online) (a) Brillouin zone of graphene. The 
circles indicate neighborhood of the K and K' points where 
the electron states have trigonal symmetry allowing for the 
photocurrent. (b) Mechanisms of photocurrent formation in 
the K valley. Solid, dashed, and dotted arrows of different 
thicknesses indicate anisotropy of velocity, optical generation, 
and scattering rate, respectively, in p-space. 



2 



of carriers in the valleys. The effective Hamiltonian de- 
scribing electron and hole states in the vicinity of the K 
and K' points has the forirP^ 
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Here, p is the momentum counted from the valley center, 

) = v v o(Px - iPy) - n(Px + Wy) 2 > ( 3 ) 

vq is the electron velocity, and /j, is the parameter of warp- 
ing that reflects the trigonal symmetry of valleys (D^h 
wave vector group). In the framework of tight-binding 
model, \x = v^a/ (A^/Sh) with a being the lattice con- 
stant P We assume that the warping is small and, there- 
fore, calculate the current to first order in /i. The energy 
spectrum of carriers in the conduction (c) and valence 
(v) bands is given by 



e ip 1=3 v oP ~ v[ip 2 cos 3</3 p , e 
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where ip p is the polar angle of the momentum p. 

Shown in Fig. lb is the energy spectrum in the K val- 
ley with the warping being included. The inequality of 
£cp' and e^2 p (as well as ei^ and £^l p ) gives rise to 
an electric current in the valley if electrons are optically 
excited from the valence to conduction band by linearly 
polarized light. In the valley K' , the warping of energy 
spectrum is opposite, Eq. Q, and the photocurrent di- 
rection is reversed. 

In the framework of kinetic theory, the photocurrent 
densities in the valleys are given by 
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where e is the electron charge, the factor 2 accounts 
for spin degeneracy, v£$, = V p e^2 are the velocities, 
and fij) and are the nonequilibrium corrections to 
the distribution functions in the conduction and valence 

and 



bands linear in the light intensity; fcp^ = /, 

fvp 1 = fv _' p due to space inversion symmetry. We con- 
sider interband optical transitions in undoped graphene 
at low temperature. Owing to electron-hole symmetry 
fcp = — /ip , and the photocurrent Eq. ^ assumes the 
formjM = 4e£ p »& ) f$. 

The steady-state correction to the distribution func- 
tion can be found from the kinetic equation 

E (KJ> & - K'l + St( ° 0) + 9% = 0, (6) 
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where W^l is the rate of intravalley electron scattering 
by static defects or impurities, the weak intervalley pro- 
cesses are neglected, gip is the optical generation rate, 
and St^ cc ^ describes electron-electron collisions. 



First, we consider the valley current in the presence 
of intensive electron scattering by impurities and neglect 
electron-electron collisions. We focus on the photocur- 
rent in the K valley and omit index v = +. In the Born 
approximation, the rate of elastic electron scattering by 
impurities W PP ' — W p > p is given by 
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where IC(q) is the Fourier component of the impurity po- 
tential correlator. The specific angular dependence of 
Wppt follows from the Hamiltonian Eq. pi). The gener- 
ation rate in the conduction band g cp is determined by 
the interband matrix elements of the velocity operator 
Vpi?p and has the form 
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Here ui is the light frequency, A/2 is the amplitude of the 
vector potential of the electromagnetic wave related to 
the intensity of incident light by I — A 2 uj 2 /(2nctfy, and 
to is the amplitude transmission coefficient, to = 2/ (n+1) 
for graphene on a substrate with the refractive index n. 

As follows from Eqs. ^ and ([6]), the valley current 
arises owing to warping-induced asymmetry in the elec- 
tron velocity v cpi the scattering rate W pp >, and the gen- 
eration rate g cp . Accordingly, to first order in fi one can 
distinguish three contributions to the current Eq. (IJ), 
X = x + x^ ^ + X^ ■ The corresponding mecha- 
nisms of the current formation are sketched in Fig. lb. 

To calculate the valley current caused by the velocity 
correction one neglects the warping in optical generation 
and scattering rates. In this mechanism, the absorption 
of linearly polarized light leads to the alignment of elec- 
tron momenta described by the second angular harmonic 
of the distribution function.^ Owing to the /^-linear cor- 
rection to the velocity, such a distribution of carriers in 
p-space implies an electric current Eq. ([!]) with 
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Here, e u = huo/2 is the kinetic energy of photoelectrons, 
T n (n = 1,2...) are the relaxation times of the nth 
angular harmonics of the distribution function, t" 1 = 
J2p> Wpp>{l — cosn8), 9 is the angle between p and p' , 
and rj = ire 2 /he is the absorbance (Ref. [3]) for normally- 
incident light. 

Another contribution to the valley current comes from 
the asymmetry of photoexcitation. Indeed, to first order 
in [i, the optical generation rate g cp contains the first 
angular harmonic, which gives rise to a photocurrent 
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The third mechanism of the current generation orig- 
inates from the asymmetry of electron scattering. The 
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solution of kinetic Eq. mn, f cp , contains the first angular 
harmonic even if the warping is neglected in the optical 
generation rate. Such a contribution to the valley current 
is given by 
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where the relaxation times are taken at the energy e u . 

Equations (|9|-(11) demonstrate that both the mag- 
nitude and excitation spectrum of the pure valley cur- 
rent are determined by the mechanisms of scattering. 
For electron scattering by unscreened Coulomb impu- 
rities in graphene, one obtains ti oc e, T2 = 3t\, and 
T3 = 5ri. Such relations yield \ ot u. In the case of scat- 
tering by short-range static defects, one has T\ oc 1/e, 
T2 = T3 = Ti/2 and, therefore, x ^ 1/w. Estimation 
shows that the valley currents in suspended graphene 
j(±) ~ 10- 4 A/cm at the light intensity 7 = 1 W/cm 2 , 
Ti = 10~ 12 s, [i = 3.6 x 10 26 g-\ and v = 10 8 cm/s. 

Now, we analyze the effect of electron-electron inter- 
action on the pure valley current. It is well known that, 
in systems with parabolic energy spectrum, the interpar- 
ticle collisions partially suppress the anisotropy of the 
distribution function. Therefore, one can expect that 
electron-electron scattering (between carriers from the 
same valley and, in particular, between carriers from dif- 
ferent valleys) can only decrease the pure valley current. 
We demonstrate below that the interparticle collisions in 
graphene may give rise to an additional contribution to 
the valley photocurrent. 

Consider the collision of a photoelectron with the mo- 
mentum p from the valley v of the conduction band with 
an electron with the momentum k from the valley v' of 
the valence band. After the collision, both electrons oc- 
cur in the conduction band with the momenta p' and fc', 
respectively, see the left inset in Fig. 2. Processes of other 
kinds are negligible in undoped graphene at low temper- 
ature and weak excitation level. The Coulomb interac- 
tion between carriers leads to the transfer of momentum 
q ~ lo/vq <C 7r/a, therefore, both electrons remain in the 
valleys they were before the collision although the valleys 
v and v 1 may be different. The momentum and energy 
conservation laws read 
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In the conic approximation, Eqs. (12 1 imply that k < p, 
p' + k' < p, and k should be antiparallel while p' 
and k' parallel to pP^ The possibility of the collisions 
is determined by corrections to the linear dispersion. 
Interaction-induced ren orma lization of the energy spec- 
trum makes it concave which forbids such Auger- 
like processes. On the contrary, the spectrum warping 
does allow for the processes. In graphene on a sub- 
strate with high susceptibility, the interaction effects are 




FIG. 2: (color online) Dependence ^(7) that determines the 
magnitude o f va lley current caused by electron-electron scat- 
tering, Eq. (17 1. Left inset sketches scattering of electrons 



with the momenta p and k into the states with momenta 
p' and k' . Right inset shows the momentum space in the 
K valley. Sectors where the Auger-like electron relaxation is 
allowed are colored, ellipses depict the alignment of photo- 
electron momenta by linearly polarized light. 



suppressed and the warping becomes prevailing. An es- 
timation shows that, for an electron with the energy 
e = 0.3 eV, this regime occurs at the effective dielectric 
constant e* > 300. Besides, interaction effects can be 
suppressed by a metallic gate or increasing the tempera- 
ture. Neglecting the interaction-induced renormalization 
of the energy spectrum, we obtain from Eqs. ( 12 1 to first 
order in p 



v (pka 2 +p , k , (3 2 ) = (13) 
-2fiv(p - p')(p' + k') \p + p' + vv\k - k')\ cos3<£ p , 



where a = ip v — ir — ipk *C 1 and f3 = ipk> — ip p i <ti 1. 
Equation (13) has solutions only for sectors of ip p where 
the second line is positive, i.e., for fj,v cos3ip p < 0. It 
means that Auger-like processes for electrons with the 
momentum p are allowed or forbidden depending on the 
sign of cos3co p , see the right inset in Fig. 2. 

The angular dependence of the electron-electron colli- 
sion rate in graphene gives rise to a pure valley current 
if the sample is illuminated by linearly polarized light. 
The mechanism of the current generation is sketched in 
the right inset in Fig. 2. The absorption of linearly po- 
larized light leads to the alignment of electron momenta 
described by the second angular harmonic of the distri- 
bution function. The photoexcited electrons are scat- 
tered by resident electrons from the valence band and lose 
their energies. Since the scattering rate is anisotropic, 
the energy relaxation leads to the formation of first har- 
monic of the distribution function and to an electric 
current if the momentum relaxation time depends on 
energy. In the case of weak electron-electron scatter- 
ing, r ee 3> Ti, where r ee is the electron-electron scat- 
tering time, the mechanism efficiency can be estimated 
as x (co) ~ ev T]T?/(euT ee ). 
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To calculate the pure valley current caused by electron- 
electron scattering we solve the kinetic Eq. ([6| with the 
linearized collision integral 
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scattering by impurities with respect to the electron- 
electron scattering rate. J(u) satisfies the boundary con- 
ditions: J(l) = J'(l) = 0, J"(l) = 2(l + 12 7 )/[ 7 (l+6 7 )], 
and J'" {I) = -36/(1 + 6 7 ) 2 . Finally, the contribution to 
the valley current induced by electron-electron scattering 
is given by Eq. ([!]) with 



ice) _ ev T}T{e u )tl 
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where u q = 2nhe 2 /(e*q) is the Fourier component of the where F ( 7 ) = j J{u) u ^du - 7 J(0)/2 + 1. The function 



Coulomb potential and the factor 4 accounts for the spin 
degeneracy. In Eq. (14), we assume that the warping 



is small and take it into account only in the arguments 
of (5-functions. The consequent simplification of St^ ee - ) 
consists in summing over the almost collinear momenta. 
Kinetic Eq. ([6| with the simplified electron-impurity col- 
lision term — /Jp /r takes the form 



9cp 




where ifk = tp P and 0(x) is the Heaviside step function. 
In the case of 1/t(e) = 2-K 2 e 4 N. l /(he* 2 e), which corre- 
sponds to the momentum relaxation time of electrons 
due to scattering by Coulomb impurities with the surface 



density Ni, Eq. ( 15 ) can be transformed into a differential 
equation 



lu {u 4 + fry)J""(u) + 7 (5u 4 - 16 7 )J"» + 
2u 2 (u 4 + 12 7 ) J'(u) + 12u(u 4 - 2 7 ) J(u) = 
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Here u = \J vop/e^, J'{u) = dJ(u)/du, and the param- 
eter 7 = irNif^vl/elj characterizes the rate of electron 



F(j) determines also the excitation spectrum of valley 
current since 7 oc 1/ui 2 . ^(7) calculated numerically 
from Eq. (16) is shown in Fig. 2. The estimation for 
huj = 1 eV and Ni — 10 12 cm -2 yields 7 ~ 10~ 2 and 
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being two orders of magnitude larger than \ 



(vel) 



X 

Thus, for graphene on a substrate with high susceptibil- 
ity, the mechanism of valley current formation caused by 
electron-electron scattering dominates. 

To summarize, we have shown that the homogeneous 
excitation of graphene with a linearly polarized light re- 
sults in a pure valley current and developed the micro- 
scopic theory of this effect. Pure valley current is accom- 
panied by no net charge current, but leads to accumula- 
tion of valley-polarized carriers at edges of the sample. 
The valley polarization breaks the time inversion sym- 
metry and also implies the local lowering of space sym- 
metry to the Z?3/j group of a single valley, which lacks 
the space inversion. The space symmetry lowering can 
be detected by optical means, e.g., by a second harmonic 
generation of the probe beam. Another possibility to 
register the valley current is to convert it into an elec- 
tric current, which can be realized in curved graphene. 
The curvature of the graphene sheet produces effective 
out-of-plane magnetic fields directed oppositely for elec- 
trons in the valleys K and K'^ Owing to the Lorentz 
force, the magnetic fields change the directions of partial 
currents in the valleys giving rise to a measurable net 
electric current. 
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